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CFMFE THERMAL ANALYSES 
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PRESSURE VESSEL AND FLUID TEMPERATURE PREDICTIONS USING 
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NODE LAYOUT FOR SINDA SPHERICAL 
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HEAT EXCHANGER, 

.0828 IN. 6061 ALUMINUM 


SIMPLIFYING ASSUMPTIONS FOR SINDA WEDGE MODELS 
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STEADY STATE TEMPERATURE 
CONTOURS FOR 22 CUBIC FOOT TANK 
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27 DEG. R. HEAT EXCHANGER 



STEADY STATE TEMPERATURE CONTOURS FOR A 200 CUBIC FOOT TANK 
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STEADY STATE TEMPERATURE 
CONTOURS FOR 200 CUBIC FOOT TANK 
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27 DEG. R. HEAT EXCHANGER 


EFFECT OF CONVECTION ON LIQUID/VAPOR INTERFACE TEMPERATURE 
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EFFECT OF CONVECTION ON LIQUID/VAPOR 

LIQUID/ 

VAPOR INTERFACE TEMPERATURE & PRESSURE 

INTERFACE 

TEMP . , 

DEG. R. 
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CONVECTION HEAT TRANSFER COEFFICIENT 
BTU/(HR-SQ. IN. - DEG. F) 



EFFECT OF WALL THICKNESS ON LIQUID/VAPOR INTERFACE TEMPERATURE 
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EFFECT OF WALL THICKNESS ON 
LIQUID/VAPOR INTERFACE TEMPERATURE 
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WALL THICKNESS, INCHES 



LIQUID/VAPOR INTERFACE TEMPERATURE VERSUS TIME 
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LIQUID/VAPOR INTERFACE 
TEMPERATURE VERSUS TIME 
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SPEAKER; GEORGE R. SMOLAK/ ANAEEX CORP. (LEWIS RESEARCH CENTER') 


John Schuster/General Dynamics Space Systems: 

The considerations on thermal mixing that you’ve done have been primarily 
limited to tanks nearly filled with liquid with a relatively small vapor ullage. 
Have you studied the opposite situation where you have mainly vapor in the tank 
and you are concerned about how that vapor will equilibrate with the liquid, and 
what the evaporation rate might be out of the liquid? 

Smolak: 

No, I didn’t go into that. The model assumes that there is a constant 
temperature at the liquid vapor interface. The model at this point does not 
include evaporation at the liquid vapor interface, that part of the 
thermodynamics is not in there. The conduction processes and convection 
processes in the vapor were included when I did include convection. That part 
of it is rigorous, but that is as far as the thermodynamics went. It is not truly 
a complete thermodynamic model by any means. 

A. A. Sonin/Massachusetts Institute of Technology: 

I wonder if you could elaborate on how you applied the convective heat transfer 
coefficient, and what the boundary conditions were on your analysis when you 
did that? 

Smolak: 

I said that I had a heat transfer coefficient of the values that were stated. The 
area over which that convective coefficient was active was the area between 
adjacent nodes in the SINDA model and the appropriate temperature difference 
was the temperature difference between the adjacent nodes. In other words, it 
was put in as a HA-Delta-T type effect. It was not done with any directionality 
or preferentially; in other words, it was done between all adjacent surfaces. If 
there is better information on how convection acts in reduced gravity, I would 
be happy to put it into the model, but I couldn’t find anything better in the 
literature. 

Robert Rudlin/Martin Marietta Denver Aerospace: 

I want to ask a question, but I observed one thing in your analysis which I 
think was very good. First, you do have to know what the convection 
coefficient or the Nusselt number is. I didn’t pick up what the Nusselt number 
was, but I did pick up that your spacial dimensions from the free surface down 
to the heat exchanger at 27 degrees, which I guess is about 3 feet, is probably 
your biggest resistor, or the primary conductor. If I draw these conclusions 
correctly, it seems to me what your saying is when you doubled everything you 
basically doubled that 3 feet to 6 feet, and, sure enough, your temperature 
difference went from 30 degrees to 60 degrees. That’s pretty much like you’d 
expect with a resistor/conductor model. Am I drawing the right conclusions in 
saying that what you’re really telling me is that you have to have your heat 
exchanger on the wall and uniformly distributed through the fluid in order to 
keep that basic dimension significantly less than 3 feet if you want to have a 
fairly well mixed fluid? In other words, don’t stir it up; just get your resistance 
path as small as possible, if you really want to control the temperature. 
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Smolak: 


That’s a long question. I would say yes, your supposition is pretty much valid. 
The scaling was done from the small tank to the larger tank and the results 
were not entirely unexpected. We did some more work that I didn’t mention 
that had to do with placing a heat exchanger on what amounted to the north 
pole in the model, which was up near the vapor bubble. That was rather 
effective in helping to destroy these temperature gradients in the fluid. Yes, 
you are on the right track. We have not explored really extensively the 
optimum positions for those heat exchangers. What we had in mind was, since 
most of the tank penetration was toward the top of the tank, that we’d 
intercept the heat coming through those penetrations with a low temperature 
fluid from the visco jet or thermodynamic vent system and that seemed to be 
rather effective. 
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